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ABSTRACT 


The  Fourier  amplitude  method  developed  in  Part  I  as  a  diagnostic  tool 
for  determining  the  sensitivity  of  the  results  of  complex  calculations  to  the 
parameters  which  enter  these  calculations  has  been  applied  to  two  chemical 
reaction  systems  involving  sets  of  coupled,  non  linear  rate  equations.  These 
were:  a)  a  5  reaction  set  describing  the  high  temperature  (6000°K)  dissociation 
of  air  and  b)  a  9  reaction  set  describing  the  constant  temperature  (2000°l<)  com¬ 
bustion  of  H2  and  0^.  We  have  evaluated  the  Fourier  amplitudes  for  all  the 
species  at  a  number  of  different  times  for  both  reaction  systems.  The  analysis 
of  these  results  verifies  the  claims  made  in  part  I.  The  relative  magnitudes 
of  the  Fourier  amplitudes  showed  a  several  order  of  magnitude  distribution  which 
permitted  a  clear  distinction  of  the  relative  sensitivity  of  the  species  concen¬ 
tration  to  uncertainties  in  the  rate  coefficients.  The  conclusions  based  on  the 
Fourier  amplitude  method  for  these  two  reaction  systems  are  in  excellent  agreement 
with  sensitivity  predictions  which  could  be  made  on  the  basis  of  previous  kinetic 
studies  of  these  systems. 


I.  INTRODUCTION 

In  the  precedes  paper,  hereafter  referred  to  as  Ip>  we  developed  a 
diagnostic  technique  to  investigate  the  sensitivity  of  the  solution  of  large 
sets  of  coupled,  non-linear  equations  to  uncertainties  or  changes  in  the  coupling 
parameters.  For  brevity,  we  will  refer  to  this  method  as  FAST,  the  Fourier 
Amplitude  Sensitivity  Test.  In  this  paper  we  wish  to  demonstrate  the  utility 
Of  FAST  by  applying  it  to  a  study  of  the  sensitivity  of  two  sets  of  coupled 
non-linear  rate  equations  describing,  respectively,  the  high  temperature 
dissociation  of  air  and  the  high  temperature  H2  -  02  reaction.  The  specific 
objective  of  the  calculations  described  here  was  to  determine  the  sensitivity 
Of  the  various  species  concentrations  at  a  particular  time  to  assumed  uncer¬ 
tainties  in  the  rate  coefficients  which  enter  into  these  calculations. 

In  Section  II  we  summarize  the  steps  which  must  be  carried  out  in  the 
application  of  FAST  to  any  system  of  coupled  non-linear  equations  and  discuss 
the  choice  of  frequencies  and  the  number  of  sampling  points  to  be  used.  In 
Section  III,  we  present  the  sensitivity  analysis  for  the  two  chemical  reaction 
systems.  Section  IV  contains  some  concluding  remarks. 

II.  GENERAL  REMARKS  ON  APPLICATION  OF  FAST 

We  wish  to  determine  the  sensitivity  of  an  “output  function",  ,n  our 
case  the  calculated  concentration  of  a  particular  species  at  a  particular  time, 
to  assumed  uncertainties  in  the  rate  coefficients  entering  the  rate  equations. 


We  write  this  output  function  as 

c(k] >  ^2’  • • • »  ^n)  =  c(£) 
where  £  is  the  n-  component  vector  of  the  r 


(2.1) 


ate  coefficients  k]t  k2,  kn> 
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The  application  of  FAST  is  carried  out  via  the  following  steps: 

Step  1)  One  must  choose  a  set  of  n  integer  frequencies,  {Wp  co2,  u  j, 

where  n  is  the  number  of  independent  rate  coefficients  which  determine  the 
value  of  the  ouiput  function. 

Step  2)  One  frequency  of  this  set  is  now  assigned  arbitrarily  to  each 
rate  coefficient  by  (see  1.2.2): 


ki 


(2.2) 
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si  n 


U)  ^  s 


(2.3) 


where  kjo)  is  the  best  estimate  of  the  rate  coefficient,  s  is  a  parameter  and 
the  u,(o'  . 

i  are  the  endpoints  of  the  estimated  ranges  of  uncertainties  of  the 
rate  coefficients  k. .  If  a  rate  coefficient  is  known  precisely,  i.e.  with  zero 
uncertainty,  then  k.  =  k!0/  which  implies  uj°)  =  0.  The  endpoints  u!°)  are 
specified  as  part  of  the  input  data  and  their  values  should  be  based  oi.  the  best 


judgement  of  the  investigator. 

Step  3)  The  output  function  c(t)  =  cdt(s)]  is  now  Fourier  analyzed  in 
s  to  obtain  the  sine  amplitudes  A^  corresponding  to  the  input  frequencies 
This  analysis  requires  the  evaluation  of  cpc(s)]  for  N  values  of  s,  where  N  is 
an  integer  which  depends  on  the  frequency  spectrum  of  c[lc(s)]. 

These  steps  will  now  be  discussed  in  more  detail. 


a)  Choice  of  Frequencies  w. 
- - - 1 


In  section  III  of  part  I  we  have  presented  a  discussion  of  the  choice 
of  the  set  of  frequencies  {w}  and  have  also  given  there  a  set  of  references 
(ref.  4)  which  contain  both  discussions  of  optimal  choices  of  sets  of  frequencies 
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as  well  as  tables  of  frequencies  which  could  be  used  for  FAST.  For  the  calcu¬ 
lations  presented  here  we  have  developed  our  own  set  of  frequencies.  We  have 
done  so  to  obtain  a  better  insight  in  how  to  handle  the  calculations  involved 
in  the  Fourier  analysis  method.  In  a  follow-on  study  we  plan  to  investigate 
whether  our  method  of  choosing  frequencies  yields  better,  equivalent  or  worse 
sets  for  use  in  FAST  than  those  presented  in  reference  4  of  part  I. 

Since  we  will  carry  out  computer  calculations,  we  are  limited  to  rational 
(or,  equivalently,  integer)  frequencies.  For  the  coupled  non-linear  equations 
under  study  here,  this  means  that  in  addition  to  the  input  frequencies  ok  and 
their  multiples,  various  linear  combinations  of  integer  multiples  of  u. 
("interference  frequencies")  will  appear  in  the  spectrum  of  c|t(s)].  This 
presents  a  problem  in  our  analysis  since  the  Fourier  sine  amplitude  of  frequency 
is  to  reflect  only  the  effect  of  the  uncertainty  in  the  value  of  the  rate 
coefficient  k£.  Clearly,  if  an  interference  frequency  coincides  with  one  of  the 
input  frequencies  ok,  say  u)£,  the  corresponding  Fourier  ampl  itude  Au  will  not 
only  reflect  changes  in  but  also  in  other  rate  coefficients.  Thus,  for 
example,  if  ^  +  co3  -  co  =  w  the  amplitude  A  is  identical  with  A,  , 

“2  (‘VVV 

and  will  therefore  reflect  not  only  changes  in  k2,  but  also  in  k]t  k3  and  k4> 

One  therefore  needs  to  choose  a  linearly  independent  set  of  input  fre¬ 
quencies  o>.  such  that 

n 

]C  Vi  *  0  ;  integer  (2.4) 

i=l 

for 

Z|ai|£M+1  (2.5) 

1=1 
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where.  M  (s  a  parameter  at  the  disposition  of  the  investigator.  He  refer  to 
such  a  set  as  being  free  of  interferences  to  order  M.  The  larger  the  chosen 
value  of  M,  the  greater  the  likelihood  that  the  Fourier  amplitude  of  each  input 
frequency  reflects  solely  the  uncertainty  of  the  corresponding  ra  e  coefficient. 
On  the  other  hand,  as  discussed  in  section  III  of  part  I  and  below,  the  larger 
the  chosen  value  of  M,  the  larger  the  maximum  value,  „max,  of  the  input  fre¬ 
quencies  which  will  still  satisfy  (2.4)  and  (2.5)  and,  correspondingly,  the 
larger  the  set  N  of  points  s  which  will  be  required  for  the  evaluation  of  the 
Founer  amplitudes.  For  tho  calculations  presented  in  this  paper  we  have  chosen 
M  =  4.  In  Table  I  we  present  sets  of  frequencies  which  are  free  of  interferences 
to  4th  order  for  systems  with  dimensionality  n  from  5  to  19,  i.e.  for  systems 
with  from  5  to  19  independent  rate  coefficients.  These  frequencies  were 
determined  via  computer  by  a  trial  and  error  procedure.  The  sets  of  frequencies 
listed  here  have  the  smallest  value  of  waax  satisfying  conditions  (2.4)  and 
(2.5)  and  are  referred  to  as  "minimal  sets", 
b)  Assignment  of  Frequencies  to  Rate  Coefficients 

The  frequences  {u^  ,  w2 . wn}  are  assigned  to  the  rate  coefficients 

{kl’  k2 . kn}  arbitrarily  since  the  results  of  these  calculations,  in  order 

to  be  useful,  must  be  independent  of  both  the  frequency  set  and  the  assignment 
of  the  frequencies.  To  check  that  the  spectrum  of  the  output  function  does  not 
contain  interfering  frequencies  which  coincide  with  the  input  frequencies,  and 
thus  affect  the  results,  one  should  reassign  the  frequencies  to  the  various  rate 
coefficients  and  repeat  the  Fourier  analysis.  If  the  calculated  Fourier  ampli¬ 
tudes  are  invariant  in  magnitude  and  sign  to  these  permutations,  then  the  set  of 
frequencies  is  free  of  interferences  to  the  order  of  the  dimensionality  of  the 
set  of  rate  equations  and  the  results  of  these  calculations  are  reliable.  If, 


on  the  other  hand,  some  of  the  Fourier  amplitudes  vary  greatly  with  the  pennu- 
tions  of  the  frequency  assignments,  a  better  set  of  input  frequencies  must  be 
generated.  By  "better"  we  refer  to  a  set  {w]  with  a  larger  parameter  M  and  a 
higher  iomax.  We  have  carried  out  such  a  permutation  analysis  for  the  -  02 
system  discussed  below  and  verified  that  our  sets  of  input  frequencies  yiel ded 
Fourier  amplitudes  ^  which  were  invariant  to  the  frequency  permutations. 

c)  Calculation  of  Fourier  Amplitudes 

As  discussed  in  section  II  of  part  I,  we  are  only  interested  in  the 
Fourier  sine  amplitudes.  For  an  input  frequency  u  ,  the  amplitude  is 

/*2it 

Aw  “  27  J  c[t(s)]  sin  u  s  ds  .  (2.6) 

To  compute  this  integral,  one  must  evaluate  the  function  c [£(s)]  at  a  set  of 
points  in  the  interval  0  ss  i  2n.  In  the  absence  of  any  information  on  the 
form  of  the  output  function  c  we  take  the  points  to  be  equally  spaced  in  that 
interval , 

s  =  *  q  =  1 ,  2,  . . . ,  N  (2.7) 


where  N  is  some  integer.  With  these  equally  spaced  points  we  can  now  conveniently 
relate  the  Fourier  coefficient  A  as  defined  by  the  integral  (2.6)  to  the  sum 

U)  _ 


A* 

w. 


N 

E 

q=l 


Sin  Sq 


c{t(sq)] 


(2.8) 


where 

A*  =  A  +  e 

0J„  (0,  (0 


(2.9) 
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The  error  term  is  given  by 


no 

EE 


'oj  "  jLu  ZmJ  A(mN-o) . ) 
a  in=l  j  J 


(2.10) 


where  Uj  is  any  frequency  in  the  spectrum  of  /[S(s)j  whi 


ch  satisfies 


mN  -  M,  =  u 
J  £ 


(2.11) 


This  error  term  is  due  to  "aliasing"  which  inevitably  occurs  whenever  a  finite 
-her  of  points  are  chosen  on  ,  interval  to  evaluate  Fourier  amplitudes.  For 

OVamnl  a  .‘r  m  ,  , 


examp, e.  if  N  equally  spaced  points  are  used,  the  amplitudes  A«>  for  . 
unavoidably  include  the  ar-nm-iwW  _  .  .  u£  1 


,  - -  •  «  i  i  Will 

unavoidably  include  the  amplitudes  of  a  component  of  frequency”*  present  in 

!  Kh’Ch  SatiSf,'e$  “*  ’■C,at''0n  —  -  -  an  arbitrary  integer, 

he  value  of  N  is  chosen  to  be  the  smallest  integer  such  that 


b.o).  f  mN 
i=l 


h i  ,m  integers 


(2.12) 


2  ^  bi‘f  ^  M  +  1 
i=l 


(2.13) 


The  Fourier  amplitudes  ^  have  been  evaluated  using  Eqs.  (2.8,  through 

If  the  output  function  c(k)  COu,d  be  expressed  as  a  polynia,  of  degree 
less  than  or  equal  to  M  in  the  variables  Uj,  i  .  , . .  where  (m#  £q. 


“j  =  In  (kpij0)) 


then  the  conditions  (2.12)  and  (2., 3)  imply  that  ,  .  o  so  that 


(2.14) 


A*  =  A 
“*  W£ 


(2.15) 


M 


The  conditions  (2.11)  through  (2.13)  also  determine  N  the  minimum  number  of  s 
points  required  for  the  evaluation  of  c[£(s)].  These  are  listed  in  Table  I. 

The  Founer  amplitudes  are  related  to  a  multiple  integral  of  an  output 
function  over  the  uncertainty  space  (see  I.  2.12).  One  way  of  obtaining  the 
Fourier  amplitudes  without  the  troublesome  interferences  and  aliasing  involved 
in  the  method  discussed  ab^/e  would  be  to  evaluate  the  n-dimensional  integral 
of  Eq.  (I.  2.12)  directly.  We  have  done  this  vi_a  simple  Monte  Carlo  integrations. 
Our  results  show  that  the  number  of  integration  points  required  by  the  Monte  Carlo 
method  to  give  a  comparable  accuracy  is  much  greater  than  the  number  required 
in  the  Fourier  analysis.  Stated  another  way,  for  a  given  number  of  integration 
points  (we  used  200),  the  degree  of  accuracy  of  the  Fourier  method  detailed 
above  far  surpasses  that  of  the  Monte  Carlo  integration, 
d)  Choice  of  the  Weight  Function 

As  emphasized  in  I,  our  final  result,  i.e.  the  relation 

A(i)  «  <dc./au0>  (2.16) 

1  * 

involves  a  weight  function  p(u.j)  which  can  be  interpreted  as  a  distribution 
function  in  u  space  which  weights  the  uncertainty  in  the  rate  coefficients.  In 
our  treatment  in  part  I,  we  used  the  weight  function  (see  Eq.  2.20) 


ai 

^ui^  "  cosh  a'. li. 


(2.17) 


l  i 


which  had  several  desirable  properties.  It  is,  for  instance,  bell  shaped  about 
u.j=0  corresponding  to  k^-k|0^,  the  "best"  value  of  the  rate  coefficient  . 

For  computational  convenience,  we  have  used  two  different  weight  functions  for 
the  calculations  presented  in  section  III  below  corresponding  to  two  different 


transformations  u,  =  ufys).  As  we  had  hoped  and  expected,  our  results  did  not 
depend  significantly  on  the  choice  of  the  weight  function.  This  indicates  that 
our  output  functions  c(k)  were  "well  behaved",  i.e.  did  not  have  any  large  local 

excursions  in  u  space.  It  is  also  in  agreement  with  the  discussion  in  the  last 
paragraph  of  section  III  of  part  I 

As  is  shown  in  part  I,  the  form  of  the  weight  function  pfu,)  depends  upon 
the  choice  of  the  transformation  function  f  of  Eq.  (1.2,3).  For  the  Fourier 
amplitude  method  to  be  at  all  useful,  our  results  on  the  relative  magnitudes  of 
the  Fourier  amplitudes  must  clearly  be  independent  of  the  explicit  form  of  the 
transformation  function.  The  insensitivity  of  our  results  to  our  choice  of 

transformation  function  (and  thus  of  the  weight  functions)  verifies  the  validity 
of  FAST  in  this  respect, 

HI.  APPLICATION  OF  FAST  TO  CHEMICAL  SYSTEMS 


We  have  applied  FAST  to  two  relatively  simple  reaction  systems  as  a 
test  of  its  utility  as  a  diagnostic  tool.  Both  systems  were  assumed,  for 
simplicity,  to  react  at  constant  temperatures  and  volume.  These  restrictions 
are  in  no  way  necessary  for  the  application  of  FAST.  The  concentrations  were 


(2) 

integrated  in  time  with  a  code  using  Gear's'  '  algorithm.  These  integrations 
were  carried  out  once  for  each  of  the  N  values  of  s  with  the  corresponding 
values  of  the  n  rate  coefficients,  t(s),  up  to  some  time  t.  This  yielded  the 
output  function  c^[lt(s),tj,  i.e.  the  concentration  of  species  i,  i=l,  2,  ...»  m, 
as  a  function  of  time.  The  Fourier  analysis  of  c^[ky(s),t]  was  then  carried  out 
at  several  selected  times.  The  uncertainties  in  the  rate  coefficients  were 
arbitrarily  assumed  to  be  much  larger  than  warranted  by  the  available  experi¬ 
mental  data  in  order  to  check  out  FAST.  It  should  therefore  be  pointed  out 
again  that  the  primary  purpose  of  the  calculations  presented  here  is  to  demonstrate 
the  validity  and  utility  of  FAST  as  a  diagnostic  tool  and  not  to  provide  new 
information  about  the  selected  reaction  systems, 
a )  High  Temperature  Air  Reaction 

We  consider  the  following  simplified  reaction  system'  ' 


02  =  m 

20  +  M 

kl*  K1 

(1) 

N2  4  M 

7 - > 

2N  +  M 

k2»  ^2 

(2) 

NO  +  M 

N  +  0  +  M 

^3*  ^3 

(3) 

n2  +  0 

NO  +  N 

k^  a 

(4) 

o2  +  n 

NO  +  0 

k5’  K5 

(5) 

The  rate  coefficients  k^ ,  i=l,  ...,  5,  refer  to  the  forward  reactions,  the  rate 
coefficients  k_^  for  the  reverse  reactions  are  obtained  from  the  equilibrium 
constant  K.  through'  the  relation  K.  =  k./k_.j.  In  the  system  and  the  one  dis¬ 
cussed  below  in  Illb,  the  equilibrium  constants  are  assumed  to  be  known  exactly, 
i.e.  with  zero  uncertainty.  The  forward  and  reverse  rate  coefficients  k^ ,  k  ^ 
for  each  reaction  therefore  have  identical  uncertainties  and  the  number  of 
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inienndent  uncertainties  is  equal  to  the  number  of  equilibrium  constants.  It 
is.  of  course,  not  necessary  to  proceed  this  way  and  one  can  instead  assign 
different  uncertainties  to  the  forward  and  reverse  rate  coefficients  which  reflect 
uncertainties  in  the  equilibrium  constants. 

The  system  was  assumed  to  react  at  a  constant  temperature  of  6000CK  and 
constant  volume.  The  rate  coefficients,  equilibrium  constants  and  initial  con- 
dmons  used  in  this  calculation  are  shown  in  Table  II. 

The  equilibrium  constants  represent  the  latest  and  “best"  literature  values  we 
could  find  and  the  rate  coefficients  listed  in  this  table  are  the  k(°)  of 
Eq.  (2.2).  In  Fig.  1  we  show  a  plot  of  the  time  evolution  of  this  reaction 
system  obtained  from  the  integration  of  the  rate  equations  using  the  rate 
coefficients  kj°)  of  Table  II. 


While  the  real  uncertainty  of  these  rate  coefficients  is  not  too  large, 
we  assumed  arbitrarily  for  the  purpose  of  testing  FAST  that  the  experimental 
uncertainties  were  i  two  orders  of  magnitude  for  each  of  the  five  independent 
sets  of  rate  coefficients.  To  take  account  of  this  uncertainty  we  then  write 


(see  Eq.  (2.2)) 

kj°)  e-4'606  <;k.  s  k!o)  e4 


606 


i  =  l. 


(3.1) 


For  our  sensitivity  analysis  we  used  the  frequency  set  (see  Table  I) 


(<4  =  (2,  42,  62,  74,  90)  (3 

for  which  the  smallest  number  N  of  evenly  spaced  points  satisfying  the  conditions 
of  Eqs.  (2.4,  2.5,  2.12,  2.13)  in  the  s  interval ,  0  <:  s  s;  2tt,  is  191. 

We  carried  out  the  Fourier  analysis  of  the  output  function  chiefs)]  at 
t  =  10  6  sec  and  10"4  sec.  From  Eqs.  (2.2),  (2.3;  and  (3.1)  we  have 

ui  =  4.606  sin  w- s 


(3.3) 


The  concentrations  of  the  various  species,  N2>  02>  NO,  0,  N,  were  calculated 
at  10  sec  and  10  sec  for  each  of  the  191  s  values  by  integrating  the  rate 
equations.  The  191  values  of  the  concentration  of  the  above  species  at  a  given 
time  define  the  output  function  ci-[t(s)]  which  is  now  Fourier  analysed  by 
computing  the  Fourier  coefficients  from  Eqs.  (2.8)  through  (2.13)  using  the  set 

of  input  frequencies  of  Eq.  (3.2).  An  example  of  the  function  c.[lt(s)]  for  NO 

-4  1 

at  t  =  10  sec  is  shown  in  Fig.  2.  This  curve  has  been  drawn  by  connecting 

the  191  calculated  concentrations  of  NO  (one  for  each  of  the  191  evenly  spaced 

values  of  s)  by  straight  lines. 

The  results  of  the  Fourier  analysis  at  t  =  10  ^  sec  are  shown  in  Table 
HI.  The  high  values  of  the  Fourier  amplitudes  associated  with  the  rate  coef¬ 
ficients  k-j ,  k_]  for  the  concentration  of  0  and  02  (see  blocks  1  and  2  of 
Table  III)  indicate  that  a  variation  in  k]  and  k_]  has  a  greater  effect  on  the 
concentration  of  0  and  02  than  the  variation  of  any  of  the  other  rate  coefficients. 
This  follows  immediately  from  Eq.  (2.16),  where  now  i  refers  to  the  0  atoms 
and  02  molecules  and  ^  =  2  refers  to  the  rate  coefficients  k]  and  k  1 .  We  can 
rephrase  this  by  saying  that,  at  this  early  time  in  the  reaction,  only  reaction 
(1)  is  important  in  the  production  of  0  atoms  and  the  disappearance  of  02 
molecules.  The  sign  of  the  amplitudes  in  the  last  column  can  be  understood  from 

Eq.  (2.16)  according  to  which  the  amplitude  is  proportional  to  (dc./du  >. 

1  £. 

Thus  the  negative  sign  of  the  Fourier  amplitude  for  k]  for  the  concentration  of 
02  (block  1,  Table  III)  simply  indicates  that  [02]  decreases  with  increasing  k] . 

The  positive  sign  for  the  Fourier  amplitude  for  k]  for  the  concentration  of  0 
atoms  (block  2,  Table  III)  indicates  that  [0]  increases  with  increasing  k-| . 

From  the  relative  magnitL-'  of  the  Fourier  amplitudes  in  blocks  %  and  5  of 
Table  III,  it  will  be  noted  that  the  formation  of  N  and  NO  at  these  early  times 
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in  the  reaction  is  most  sensitive  to  reactions  (1)  and  (4).  This  is  quite 
reasonable  since  at  t  =  ICf6  sec  the  formation  of  N  and  NO  via  reactions  (2) 
and  (5)  cannot  play  any  important  role  owing  to  the  small  value  of  k(o)  as 
compared  to  k<°>  (see  Table  II).  A  number  of  other  conclusions  as  to  the 
sensitivity  of  the  concentrations  of  the  various  species  with  respect  to  vari¬ 
ation  of  the  rate  coefficients  can  be  drawn  from  a  study  of  the  results  in  Table 

III.  All  these  conclusions  are  in  accord  with  one's  knowledge  and  intuition  of 
this  rather  simple  reaction  system. 

In  Table  IV  we  present  our  results  for  t  »  10'4  sec.  At  this  time,  as 
can  be  seen  from  Fig.  1,  the  concentrations  of  reactants  and  products  are, 

except  for  Nj,  significantly  different  from  their  initial  values.  The  coupling 

between  the  various  reactions  has  come  in+n  ni, 

lnt°  Play  “inch  more  strongly  at  this  later 

time  and  one  would  expect  that  the  influence  of  the  uncertainties  in  the  rate 
coefficients  on  the  concentrations  of  the  various  species  is  more  complex.  This 
can  readily  be  verified  from  the  data  in  Table  IV,  where  there  is  now  no  longer 
such  a  pronounced  order  of  magnitude  effect  in  the  Fourier  amplitudes.  Let  us 
arbitrarily  use  a  factor  of  about  10  to  distinguish  between  the  "relatively 
important"  amplitudes  (i.e.  the  concentration  of  the  species  is  very  sensitive 
to  the  change  of  a  specific  rate  coefficient)  and  the  "relatively  unimportant" 
amplitudes  (i.e.  the  concentration  of  the  species  is  not  very  sensitive  to  the 
change  of  a  specific  rate  coefficient).  Then  the  data  of  Table  IV  indicates 
that  [02]  and  [0]  are  affected  primarily  by  uncertainties  in  k, ,  [N  ]  is 
affected  about  equally  by  uncertainties  in  all  of  the  rate  coefficients,  [N]  is 
least  affected  by  uncertainty  in  k, ,  and  [NO]  is  least  affected  by  uncertainty 
in  k2.  The  least  "sensitive"  rate  coefficient  is  k2>  which  determines  the  rate 
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is  smaller  than  the 


of  dissociation  of  .  This  is  not  surprising  since  k,/0^ 
other  rate  coefficients  hy  2  to  4  orders  of  magnitude.  The  important  point  to 
note  is  that  owing  to  the  strong  coupling  of  all  the  reactions  at  this  stage  of 
the  kinetic  development,  a  number  of  the  elementary  rate  processes  contribute 
significantly  to  the  production  and  removal  of  most  of  the  species,  so  that  the 
concentration  of  these  species  are  quite  sensitive  to  uncertainties  in  several 
rather  than  just  one  rate  coefficient.  This  result  is  certainly  not  unexpected. 

b)  The  H2  -  02  Reaction 

The  H2  -  02  combustion  system  was  modeled  by  the  following  set  of  chemical 
reactions^ : 


H2  +  M 

=*  2H  +  M 

kl ,  K1 

(1) 

0 2  +  M  P~ 

-  *  20  +  M 

k2 ,  ^2 

(2) 

OH  +  M 

*  0  +  H  +  M 

k3,  K3 

(3) 

H20  +  M  ^ 

H  +  OH  +  M 

k4»  K4 

(4) 

H  +  02<F= 

—  OH  +  0 

k5’  K5 

(5) 

0  +  H2£= 

=±  OH  +  H 

k6  ’  K6 

(6) 

H2  +  OH  «= 

=*  H20  +  H 

k7  ’  K7 

(7) 

o2  +  h2  «= 

—  20H 

k8»  Kg 

(8) 

OH  +  OH  f= 

=*  H20  +  0 

kg  ,  Kg 

(9) 

The  system  was  assumed  to  react  at  a  constant  temperature  of  2000°K  and  constant 
volume,  with  only  H^  and  02  present  at  t  =  0. 

The  initial  conditions,  rate  coefficients  and  equilibrium  constants  are 
shown  in  Table  V.  The  rate  coefficients  listed  in  this  table  are  the  k|0^  of 
Eq.  (2.2);  the  equilibrium  constants  are  assumed  to  be  known  with  zero  uncertainty. 


We  assumed  arbitrarily  that  the  experimental  uncertainties  for  each  of 


the  nine  independent  rate  coefficients  were  *  one  order  of  magnitude  so  that 


-2.303 

e 


£k(?)  e2*303 


i  *1,  2 .  9 


(3.4) 


For  our  analysis  we  used  the  frequency  set 


M  =  (19,  59,  91  ,  113,  133,  143,  149,  157,  161  )  (3.5) 

for  which  the  smallest  number  N  of  evenly  spaced  points  in  the  s  interval, 

0  £  s  £  2rr ,  is  630. 

A  plot  of  the  concentrations  of  the  various  species  for  the  time  interval 
10"4  £  t  s:  10  seconds,  calculated  with  the  kj°)  of  Table  V  is  shown  in  Fig.  3. 

We  have  employed  FAST  to  study  the  sensitivity  of  the  species  concentrations  to 
uncertainties  in  the  rate  coefficients  at  10~8  seconds  and  10"3  seconds.  These 
results  are  shown  in  Tables  VI  and  VII. 

At  the  very  early  time  of  10"8  seconds  (Table  VI),  the  H2  and  02  concen¬ 
trations  have  not  changed  significantly  from  their  original  value  (to  within 
8  significant  places)  and  the  Fourier  amplitudes  are  merely  round-off  error. 

For  H  atoms,  reaction  (1)  seems  to  be  the  major  contributor  to  its  formation 
at  this  early  time  and  its  rate  coefficient  is  the  most  sensitive  in  determining 
the  accuracy  of  [H].  Reactions  (7)  and  (8),  which  together  also  produce  H  atoms, 
are  also  important  but  the  concentration  of  H  atoms  is  not  as  sensitive  to  ky 
and  k8  as  it  is  to  ky  The  sensitivity  of  [H]  to  these  two  rate  coefficients  has 
identical  values  which  is  not  surprising  since  [H]  is  proportional  to  k7-kQ  at 

this  early  time.  The  same  argument  also  holds  for  the  production  of  H20  since 

—0 

its  concentration  at  10  seconds  is  proportional  to  k^kg.  The  H20  concentration 
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IS  also  sensitive  to  kg,  since  reaction  (5)  appears  to  be  partly  responsible 
tor  the  production  of  the  OH  which  enters  into  reaction  (7).  The  concentration 
of  0  atoms  is  most  sensitive  to  k2  which  is  again  as  expected.  Finally,  the  OH 
concentration  is  most  sensitive  to  kg  which  is  also  most  reasonable,  since  at 
this  early  time,  only  the  direct  reaction  between  02  and  H2  would  be  expected 
to  effect  the  OH  concentration.  Again,  the  results  of  FAST  are  in  good  agreement 
with  what  could  be  predicted  on  the  basis  of  chemical  knowledge. 

At  t  =  10  seconds  (Table  VII)  the  story  is  quite  different.  As  can  be 
seen  from  Fig.  3,  significant  changes  from  the  initial  concentrations  have 
occured  by  this  time,  and  the  reactions  are  now  strongly  coupled.  This  is 
reflected  in  the  array  of  Fourier  amplitude  in  Table  VII  as  compared  to  those 
of  Table  VI.  Now,  for  instance,  rate  coefficient  kg  corresponding  to  the  chain 
branching  step  (5)  which  produces  H  and  0  atoms  is  the  most  influential  one  in 
controlling  the  concentration  of  all  species.  The  H  atoms  concentration,  in 
addition,  is  equally  sensitive  to  k7>  All  other  rate  coefficients,  except  for 
k3,  which  is  very  small  for  all  species,  have  about  equal  sensitivities  to 
within  a  factor  of  10  as  measured  by  their  Fourier  amplitudes.  According  to 
our  theory  and  calculations  then,  the  rate  coefficient  which  needs  to  be 
determined  with  the  highest  accuracy  for  the  accurate  calculation  of  all  species 
concentration  at  10  3  seconds  is  kg  with  k^  next  in  line.  Uncertainties  in  the 
other  rate  coefficients  will  have  a  smaller  effect  on  the  specier  concentrations. 


iv.  CONCLUDING  REMARKS 

It  should  be  realized  that  the  applications  presented  here  are  only  a 

first  cut  at  the  subject.  It  appears  that  FAST  is  a  useful  diagnostic  tool  in 

determining  the  sensitivity  of  the  results  of  complex  calculations  to  the 

parameters  which  enter  the  calculations.  The  Fourier  amplitude  method  clearly 

needs  further  testing  on  more  complicated  systems  than  the  ones  dealt  with  in 

this  paper.  Such  tests  should  also  involve  comparison  with  various  types  of 
brute  force"  calculations. 

In  addition  to  applications  to  larger  systems,  there  are  also  a  number  o 

theoretical -computational  questions  which  need  further  investigation  Some 

examples  arei  the  influence  of  the  range  of 'uncertainties,  with  different 

uncertainties  for  different  coupling  parameters,  on  the  relative  magnitudes  of 

the  Fourier  amp,itudes;  the  influence  of  the  choice  of  the  input  frequencies  ,„> 

aid  the  spacing  of  the  N  values  of  the  parameter  s  on  the  Fourier  amplitudes;  an 

the  question  whether  the  absolute  magnitude  of  the  Fourier  amplitudes  can  be 

used  as  a  predictor  for  the  actual  change  of  the  output  function  for  a  given 

Change  of  the  coupling  parameter.  Additional  work  is  being  carried  out  to  answe, 
these  questions. 
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FIGURE  CAPTIONS 


Fig.  1.  Time  evolution  of  the  N2  -  02  system  based  on  the  kinetic  data 
in  Table  II. 

Fig.  2.  The  function  c.[£(s)]  for  NO  at  10"^  sec  between  0  and  2rr  as 
a  function  of  s.  This  curve  has  been  drawn  by  connecting  the 
191  calculated  values  of  the  concentration  of  NO  (corresponding 
to  the  191  values  of  s)  by  straight  lines. 

-4 

Time  evolution  of  the  H2  -  02  systems  between  t  =  10  sec  to 
t  =  iO  sec  based  on  the  kinetic  data  in  Table  V. 


Fig.  3. 
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Table  I 


FREQUENCY  SETS  AND  CORRESPONDING  NUMBER  OF 
POINTS  AVOIDING  INTERFERENCES 
THROUGH  4th  ORDER 


Dimension  (n) 

Frequency  Set  {u} 

Minimum  Number  of  Points  (N ) 

5 

11  ,  21  ,  27,  35,  39 

142 

5* 

2,  42,  62,  74,  90 

191 

6 

1 ,  21 ,  31,  37,  45,  49 

182 

6* 

2,  42,  62,  74,  90,  98 

231 

7 

17,  39,  59,  69,  75,  83,  87 

334 

8 

23,  55,  77,  97,  10/,  113, 

121,  125 

486 

9 

19,  59,  91  ,  113,  133,  143 

149,  1  57,  161 

630 

10 

25,  63,  103,  135,  157,  177, 
187,  193,  201,  205 

806 

11 

41  ,  67,  105,  145,  177,  199, 
219,  229,  235,  243,  247 

974 

12 

31  ,  87,  113,  151  ,  191  ,  223, 
245,  265,  275,  281,  289,  293 

1158 

13 

23,  85,  141  ,  167,  205,  245 
277,  299,  319,  329,  335,  343, 

347 

1374 

14 

87,  133,  195,  251,  277,  315, 
355,  387,  409,  429,  439,  445, 
453,  457. 

1814 

15 

67,  143,  189,  251,  307,  333, 
371,  411,  443,  465,  485,  495, 
501  ,  509,  513 

2038 

16 

73,  169,  245,  291,  353,  409, 
435,  473,  513,  545,  567,  587, 
597,  603,  611,  615 

2446 

17 

85,  145,  241  ,  317,  363.  425, 

'  481  ,  507,  545,  585,  617,  639, 
659,  669,  675,  683,  687 

2734 

18 

143,  229,  289,  385,  461,  507, 
569,  625,  651,  689,  729,  761, 
783,  803,  813,  819,  827,  831 

3310 

19 

149,  275,  361,  421,  517,  593, 
639,  701  ,  757,  783,  821,  861, 
893,  915,  935,  945,  951,  959, 

963 

3848 

*Thcse  sets  of 

5  and  6  frequencies  are  not  minimal, 

but 

have  been  used  and  are 

only  included  for  completeness.  They  are  also  free  of  interferences  to  order 


less  than  5, 


- . — 


23 


Table  II 


INITIAL  CONDITIONS,  RATE  COEFFICIENTS  AMD  EQUILIBRIUM 
CONSTANTS  FOR  THE  HIGH  TEMPERATURE  AIR  REACTION 

Initial  Conditions 
LN2]  -  8  x  10-8  moles/cc 
[02]  =  2  x  10'8 


Rate  Coefficients 


k{°)  -  8.5  x I O1 0 

(moles/cc)"^ 

k(°)  =  3.0  x  !07 

II 

=  8.0  x  109 

II 

k|0^  =  9.0  x  1010 

II 

kg0^  =  8.0  x  1011 

II 

Equilibrium  Constants 

K-j  =  7.8  x  10~^  (moles/cc) 

K2  =  1.0  x  10‘7 

II 

K3  =  1.3  x  10-5 

H 

K4  =  9.0  x  10"3 

- 

K,  =  5.9  x  101 

1 


M-faiin .I- 1,  if 


Table  III 


FOURIER  AMPLITUDES  OF  THE  CONCENTRATIONS  AT  10"6  SEC 
FOR  THE  HIGH  TEMPERATURE  AIR  SYSTEM 


Input  Frequency 


Rate  Coefficient 


Fourier  Amplitude 


[02]  x  109  moles/cc 


-4.47 


-0.0244 


-0.0227 


-4.7  x  10"4 
-3.9  x  10"4 


[0]  x  109  moles/cc 


-0.36 


0.018 


-0.029 


0.0027 


[^2^3  x  109  moles/cc 


7 


* 


Table  IV 


A 

FOURIER  AMPLITUDES  OF  THE  CONCENTRATIONS  AT  10“4  ^EC 
FOR  THE  HIGH  TEMPERATURE  AIR  SYSTEM 


Input  Frequency 

Rate 

Coefficient 

Fourier  Amplitude 

X 

r— i 
Csj 
O 
i _ i 

108  moles/cc 

2 

kl 

-1.1 

90 

k5 

-0.147 

74 

k4 

-0.087 

62 

k3 

-0.069 

42 

k2 

-0.044 

[0]  x  108  moles/cc 

2 

kl 

2.17 

.  62  - 

k3 

0.228 

90 

k5 

0.164 

42 

k2 

0.078 

74 

k4 

-0.00146 

[Nz]  x 

g 

10  moles/cc 

74 

k4 

-2.6 

2 

kl 

-2.3 

62 

k3 

-1.0 

90 

k5 

-0.70 

42 

k2 

-0.55 

[N]  x  10^  moles/cc 


74 


3.3 


Table  V 


Initial  Conditions,  the  Rate  Coefficients  and  Equilibrium 
Constants  for  the  H p  -  Op  Reaction 


Initial  Conditions 

[Up]  =  8.0  x  10"^  moles/cc 
[Op]  -  4.0  x  10"^  moles/cc 


Rate  Coefficients 

k°  =  5.78  x  IQ4  (moles/cc)"1 sec 

k°  •=  4.47  x  103 

k°  =  1 .03  x  10b 

k°  -  6.79  x  105 

k°  =  1  .23  x  1012 
Equilibrium  Constants 

K]  =  1.59  x 
Kp  =  2.69  x 
K3  =  1.14  x 
=  1 .63  x 
K5  =  2.37  x 


10"11  moles/cc 
-12  „ 


10 


10 


•11 


10 


-12 


10 


-1 


kg  =  4.0  x  1013  (moles/cc)_1sec 
k^  -  1.51  x  1013 
k°  -  2.20  x  106 
kg  =  1.6  x  1014 

K6=!.4  - 

K?  =  9.76  - 

Kg  =  2.02  - 

Kg  =  6.76  - 


21 


Table  VI 


FOUR  1 CR  AMPLITUDES  OF  THE 
FOR  THE  IL 


CONCENTRATIONS  AT  10"8  SEC 
-  0?  REACT] OR 


Reproduced  from 
best  f.vailable  copy. 


Rate  Coefficient 


Fourier  Amplitude _ Rato  Coefficient _ Fourier  Ampl 


[h2] 

Insensitive  to  uncertainties 
in  any  rate  constant. 

[H]  x  1019 


[OH]  x  10 


18 


[0  2] 

Insensitive  to  uncertainties 
in  any  rate  constant. 


[0]  x  10 


20 


kl 

4.67 

k2 

1.79 

k8 

0.202 

k6 

-0.016 

k7 

0.202 

kl 

0.0065 

k5 

0.010 

k5 

0.0055 

k6 

0.0023 

k3 

2.25x1 0"5 

k2 

0.00166 

k4 

-1 .95x1 O-6 

k4 

8.0x1 0-5 

k9 

7.87x1 0-6 

k9 

4.1x1 0-5 

k8 

1 . 20x1 0~6 

k3 

-3.0xl0~5 

k7 

-8.46x1 O-7 

[H20]  x  10 
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k8 

5.9 

k8 

2.03 

k7 

-0.0203 

k7 

2.03 

k9 

0.00176 

k5 

0,121 

k5 

•0.0014 

k6 

0.0066 

k2 

1.61xl0"4 

k4 

8.35x1 0“4 

k6 

9.73x1 0“5 

kl 

-6.35x1 0"4 

kl 

7.08x1 0"5 

k9 

4.22x1 0-4 

k4 

-7.43xl0"6 

k3 

-2.81 xl 0"4 

k3 

3. 13x1  O''6 

k2 

1 .33x1 0"4 

Table  VII 


FOR  THE  H2  -  02  REACTION 


Rate  Coefficient  Fourier  Amplitude 

Rate  Coefficient  Fourier  Amplitude 

[H2]  x  109 

[H]  x  109 

k5 

-3.71 

k5 

1.92 

k7 

-0.346 

k7 

1.92 

k6 

-0.313 

k6 

0.122 

k8 

-0.172 

k8 

0.102 

k2 

-0.108 

k2 

0.061 

k4 

-0.101 

k9 

0.0466 

k9 

-0.0905 

k4 

0.0439 

kl 

-0.0491 

kl 

0.0176 

k3 

0.0119 

k3 

-0.00318 

[02J  x  109 

i — i 

o 

i _ i 

X 

o 

o 

k5 

-1.77 

k5 

5.77 

k7 

-0.153 

k7 

0.420 

k6 

-0.089 

ki 

0.164 

k8 

-0.070 

k6 

0.162 

k2 

-0.047 

k9 

0.145 

k4 

-0.047 

k2 

0.116 

k9 

-0.044 

k8 

0.115 

k‘I 

-0.031 

k4 

0.10 

k3 

0.0091 

k3 

-0.061 

[OH]  x  1010 

[H20]  x  109 

k5 

4.63 

k5 

2.52 

Ir  n  001  i. 
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